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The PHENIX experiment has measured the production of tt°s in Au+Au collisions at -^s NN = 
200 GeV. The new data offer a fourfold increase in recorded luminosity, providing higher precision 
and a larger reach in transverse momentum, px, to 20 GeV/c. The production ratio of rj/n a is 
0.46±0.01(stat)±0.05(syst), constant with pr and collision centrality. The observed ratio is consis- 
tent with earlier measurements, as well as with the p+p and d+Au values. ir° are suppressed by 
a factor of 5, as in earlier findings. However, with the improved statistical precision a small but 
significant rise of the nuclear modification factor Raa vs pr, with a slope of 0.0106±q;oo29 [Gev/c] _1 , 
is discernible in central collisions. A phenomenological extraction of the average fractional parton 
energy loss shows a decrease with increasing pr- To study the path length dependence of sup- 
pression, the 7T° yield was measured at different angles with respect to the event plane; a strong 
azimuthal dependence of the n° Raa is observed. The data are compared to theoretical models of 
parton energy loss as a function of the path length, L, in the medium. Models based on pQCD are 
insufficient to describe the data, while a hybrid model utilizing pQCD for the hard interactions and 
AdS/CFT for the soft interactions is consistent with the data. 

PACS numbers: 25.75. Dw 



I. INTRODUCTION 



Discovery of the suppression of high transverse momen- 



tum (pr) hadrons in relativistic heavy ion collisions [1|-3| 
and the absence of such suppression in dAu collisions 4| 
inspired intense theoretical work during the past decade. 
The phenomenon was immediately interpreted, in fact, 
even predicted 5- 7], as the energy loss of a hard scattered 
parton in the hot, dense strongly-interacting quark- gluon 
plasma (QGP) formed in the collision. Prompted by the 
large amount of very diverse experimental data from the 
Relativistic Heavy Ion Collider (RHIC) - namely, by sup- 
pression patterns at various collision energies, colliding 
systems, and centralitics - several models have been de- 
veloped, based mostly on perturbative quantum chromo- 
dynamics (pQCD) (see Section UlI El as well as [8|). The 
suppression patterns are quantified by the nuclear mod- 
ification factor -Raa, defined for single-inclusive 7r°s as 



Raa{pt) 



(l/N^)d 2 NlJd PT dy 
(T A b> x d 2 af p /dp T dy ' 



(1) 



where cr^ ° is the production cross section of n° in p+p col- 
lisions, (Tab) = (Aran) /fJpp 1 is the nuclear overlap func- 
tion averaged over the relevant range of impact parame- 
ters, and (N co \i) is the number of binary nucleon-nucleon 
collisions computed with <r " el . Despite their different ap- 
proaches, several models |9Hl2| were able to describe the 
Pt and centrality dependence of Raa within experimen- 
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tal uncertainties. At the same time, those models pro- 
vided very different estimates of medium properties such 
as the transport coefficient q, the average 4-momentum- 
transfer-squared per mean free path of the outgoing par- 
ton within the medium. For this reason, Raa alone does 
not provide sufficient constraint for extracting medium 
properties such as q from the theoretical predictions, be- 
cause it averages the varying energy losses along many 
different paths of the parton in the medium. 

While dihadron correlation measurements are a suc- 
cessful approach to constrain (L) of the parton in the 
medium [13J , the single particle observable Raa typically 
has smaller statistical errors and higher px reach. In ad- 
dition, if Raa is measured as a function of the azimuthal 
angle with respect to the event plane of the collision, the 
average path length (L) can be constrained [13, EH ■ In 
all but the most central ion-ion collisions, the overlap 
region of the nuclei is not azimuthally isotropic. The av- 
erage distance the parton traverses before emerging and 
fragmenting varies as a function of the angle with respect 
to the event plane. Each collision centrality A<fi class se- 
lects different (L) values, so the differential observable 
-Raa(A^) directly probes the path length dependence of 
the energy loss. 

The first measurements of azimuthal as ymm etries of 
nuclear suppression and collective flow [14j - |l6j ] used 7r°s 
as the probe, which has the advantage that 7r°s are rel- 
atively easy to identify over a very wide px range in a 
single detector - a crucial factor in mitigating systematic 
uncertainties. As pointed out in jlEj . both collective flow 
and azimuthal dependence of nuclear suppression can be 
formally defined at any px, but they have historically and 
conceptually different roots. The notion of collective flow 
originates in lower px phenomena, and is usually inter- 
preted as a boost to the original px spectrum (of partons 
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or final state particles) in the direction of the highest 
pressure gradient. In contrast, Raa and -Raa(A</>) are 
typically used to describe high px behavior, and their 
decrease from unity interpreted as a loss of parton mo- 
mentum due to the presence of a medium. In this paper, 
results on tt° production, the nuclear modification fac- 
tor i?AA, and its azimuthal dependence in terms of the 
event-plane-dependent -Raa(A0) are presented. The re- 
sults presented here are based upon the data collected in 
the 2007 RHIC run. The data sample is four times larger 
than that of [HI, [l7| ■ The dedicated reaction plane de- 
tector installed in 2007 offers improved event-plane 
resolution. 



II. EXPERIMENTAL DETAILS 



collisions iV co n, and impact parameter b are estimated 
(see Table H|). 



TABLE I: Average iVpart , N co \\, impact parameter and par- 
ticipant eccentricity [221 ] for all centrality classes. 



Ccntrality 


(JVpart ) 


(N coU ) 


W 


(epart) 


(%) 






[fm] 




00-10 


325.8±3.8 


960.2±96.1 


3.1±0.1 


0.105 ± 0.004 


10-20 


236.1±5.5 


609.5±59.8 


5.6±0.2 


0.198 ± 0.008 


20-30 


167.6±5.8 


377.6±36.4 


7.3±0.3 


0.284 ± 0.010 


30-40 


115.5±5.8 


223.9±23.2 


8.7±0.3 


0.358 ± 0.011 


40-50 


76.2±5.5 


124.6±14.9 


9.9±0.4 


0.425 ± 0.013 


50-60 


47.1±4.7 


63.9±9.4 


10.9±0.4 


0.495 ± 0.016 


60-70 


26.7±3.7 


29.8±5.4 


11.9±0.5 


0.575 ± 0.023 


70-80 


13.7±2.5 


12.6±2.8 


12.6±0.8 


0.671 ± 0.024 


80-93 


5.6±0.8 


4.2±0.8 


13.9±0.5 


0.736 ± 0.021 



A. Data set 



B. Reaction plane 
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FIG. 1: (Color online) PHENIX experimental setup in the 
2007 data taking period. 

This analysis used 3.8 x 10 9 minimum bias Au+Au col- 
lisions at ^/sjvjv = 200 GeV recorded by the PHENIX ex- 
periment [19| at RHIC in 2007. The experimental setup 
is shown in Figure [TJ Collision centrality was determined 
from the amount of charge deposited in the Beam-Beam 
Counters (BBC, 3.0 < \rj\ < 3.9). From a Monte Carlo 
calculation based on the Glauber model [I3,[2lj], the aver- 
age number of participants -/Vp art , the number of binary 



Each noncentral nucleus-nucleus collision has a well- 
defined reaction plane, given by the beam direction and 
the impact parameter vector of the actual collision. Al- 
though this reaction plane cannot be directly observed, 
an event plane can be experimentally determined event- 
by-event using the method discussed in detail in [23T ]. 

In order to reduce the biases to the event plane de- 
termination from physical correlations such as Hanbury- 
Brown-Twiss (HBT), resonance decay, and especially 
high-pT jet production, it is necessary that the event 
plane is determined with a large rj gap with respect to 
the high-pT measurement 24|. Therefore, in this anal- 



ysis measurements from two detectors were combined, 
located along the beam direction to the North and South 
of the interaction region. The first is a pair of Muon 
Piston Calorimeters (MPC) [2j| covering 3.1 < \rj\ < 3.9 
in pseudorapidity and consisting of 240 2.2 x 2.2 x 18 cm 3 
PbW04 crystals each. The second is a pair of Reaction 
Plane Detector (RxNP) [l8j |. which are plastic scintilla- 
tors, with 20 mm of lead converter in front of it. The 
RxNP is divided into 12 azimuthal segments and further 
divided radially into outer (RxNPout) and inner (RxN- 
Pin) rings. The outer ring covers 1.0 < |?7| < 1.5 and the 
inner ring covers 1.5 < |?7| < 2.8. The current analysis 
used RxNPin only. The method to establish the event 
plane from the combined MPC-RxNPin information is 
identical to the one used in . 

In order to estimate the resolution of the event plane, 
it is measured independently by the north and south 
detectors, and ^>St respectively. The resolution is 
then characterized by (cos[2(*Ar — ^s)])- Higher values 
indicate better resolution. The resolution is centrality- 
dependent, as shown in Figure [5J 



C. Neutral pions 

Neutral pions are measured via the ir° — > 77 decay 
channel. Photons are identified in the PHENIX Electro- 
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FIG. 2: (Color online) Event plane resolution as a function of 
collision centrality expressed in terms of N pa ,rt , using only the 
BBC, and using the combined MPC and RxNPin detectors. 
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FIG. 3: (Color online) Invariant mass spectrum of two pho- 
tons (black) and the corresponding mixed events (red) at 
7 < Pt < 7.5GeV/c in minimum bias collisions. Vertical 
lines indicate a ±2.5 a integration window. 



magnetic Calorimeter (EMCal, described in [2g|) consist- 
ing of two subdetectors, both extending to \r\\ < 0.35 in 
pseudorapidity and are located at 5.1m radial distance 
from the collision point. The analysis uses data from 
the lead-scintillator (PbSc) sampling calorimeter, which 
comprises six sectors covering 3/8 of the full azimuth 
and has a 5.5x5.5 cm 2 granularity and depth of 18 radi- 
ation lengths. Photons are identified using various cuts 
on the shower shape observed in the calorimeter as well 
as by comparing the observed shapes to an ideal one, 
parametrized using well-controlled test beam data [26| . 
Since this analysis is restricted to the pt region above 



5GcV/c, the hadron contamination is small; hadrons in 
this energy region typically deposit only a small fraction 
of their energy in the EMCal. 



The invariant mass m 77 is calculated in bins of pho- 
ton pair pt from each pair of photons, provided the 
pair passes the energy asymmetry cut a < 0.8 where 
a = \E. 
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Ey 2 \/(Ey 1 + E l2 ), and the distance between 
the impact positions of the two photons is larger than 
8 cm. An example m 77 distribution is shown in Fig. [3] 
For the event-plane-dependent studies the procedure is 
repeated in six 15°-wide bins of angles Ac/> with respect 
to the event plane. The combinatorial background is 
estimated with the event mixing technique where pho- 
tons from one event are combined with photons from 
other events, which satisfy the same global conditions 
(vertex position, centrality, event plane direction), and 
calculated. The mixed-event distributions 
are then normalized and subtracted from the real event 
distributions. The resulting 7r° peaks are a =10-11 MeV 
wide, depending on centrality, and have very small resid- 
ual background due to the inherent correlations in real 
events not reproducible by the mixed-event technique. 
This residual background is fitted to a second-order poly- 
nomial in the regions below and above the tt° peak. This 
polynomial shape is then subtracted from the m 77 dis- 
tribution. The raw ir° yields are extracted by integrating 
the resulting histogram in a ±2.5<r wide window. 



In order to establish the combined effects of accep- 
tance and 7T° detection efficiency, single 7r°s are gener- 
ated with a distribution uniform in 4> & n d extending to 
1 77 1 < 0.5 in pseudorapidity, then simulated in the full 
GEANT3 US framework of PHENIX. After the GEANT3 



output is tuned to reproduce the inactive detector areas 
as well as the peak positions and widths observed in real 
data, the simulated 7T°s are embedded into real events. 
The embedded output can then be analyzed with the very 
same tools as the real events. 



At high pt, the two decay photons may be so close 
that the EMCal can no longer resolve them as two parti- 
cles and provide the proper energies and impact points. 
The two photons "merge" into one cluster, and the cor- 
responding 7T° cannot be reconstructed from m 77 . Such 
merged clusters were rejected by various shower profile 
cuts, and the loss was determined by simulated 7r°s em- 
bedded into real events and analyzed with the same cuts. 
At 11 GeV/c merging happens only for the most symmet- 
ric decays resulting in a 5% loss of 7r°s. At 17GeV/c the 
correction is 50%. At p T = 20GeV/c about 70% of tt°s 
are lost due to this effect. The systematic uncertainties 
were estimated by comparing ir° yields extracted in bins 
of asymmetry (a). The ir° yields are corrected for the pt 
bin width by fitting the invariant yield to a power-law fit 
and adjusting the yield to correspond to the one at the 
center of the pt bin. 
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TABLE II: Typical (minimum bias) values of systematic un- 



certainties of the 


invariant yields of ir°. 








„ rr^^A/ I,, l 

pT^CV/C J 


indep. 


6 


8 


1U 


ID 


lype 


i leia extr. \/o) 




5.0 


4.0 


6 . U 


Z.U 


r> 
JD 


E scale (%) 




6.0 


6.0 


7.0 


7.0 


B 


PID (%) 




4.0 


3.0 


4.0 


5.0 


B 


Merging (%) 








4.5 


28.0 


B 


Acceptance (%) 


1.0 










B 


Off-vertex (%) 


1.5 










C 


Total (%) 


1.8 


8.8 


7.8 


9.7 


29.4 





D. Systematic uncertainties 

Systematic uncertainties are characterized as follows. 
Type A uncertainties are point-to-point uncorrelated 
with pt- Type B uncertainties have point-to-point cor- 
relations that cannot be characterized by a simple multi- 
plicative factor, but vary smoothly with px- Finally, type 
C uncertainties would move all points up or down by a 
common multiplicative factor, a typical example being 
the uncertainty on N co \\ in Raa- 

The type B systematic uncertainty of the tt° raw yield 
extraction has been estimated by comparing yields ob- 
tained in windows of varying widths. The uncertainty 
is less than 5% for peripheral collisions (low multiplicity, 
small combinatorics) and reaches about 7% in central 
collisions. 

The uncertainty on the efficiency of the photon identi- 
fication (PID) is estimated comparing fully corrected tt° 
yields obtained with various PID cuts. The uncertainty 
is 2-4% at 5-8GeV/c, and increases both with centrality 
and with px- It is of type B. 

The uncertainty on the energy scale is estimated from 
how well the peaks and widths of simulated ir°s embed- 
ded in real events agree with the measured peaks and 
widths at each centrality. The difference is less than 1% 
at 5-8 GeV/c. Due to the steeply falling 7r° spectrum this 
less-than-1% uncertainty of the energy scale translates to 
about 7% uncertainty on the ir° invariant yield. 

The uncertainty due to the photon-merging correction 
is estimated as follows. Raw yields at high px are ex- 
tracted in different asymmetry windows both from real 
data and simulated decay photon pairs embedded in real 
data. Apart from small and precisely calculable accep- 
tance effects, the true asymmetry distribution is flat, and 
at any given px one should observe the same raw ir° 
yield, for instance, in the window 0.4 < a < 0.6 and 
0.6 < a < 0.8. However, lower asymmetry means a 
smaller opening angle of the decay leading to a greater 
probability for the photons to merge. Therefore, the mea- 
sured asymmetry distribution at high px is not flat. To 
determine the photon-merging correction and its system- 
atic uncertainty, a series of raw yield ratios in different 
asymmetry bins were compared between data and sim- 
ulation. The uncertainties on the ir° spectra due to the 
merging correction are px and centrality dependent. 



The uncertainty due to acceptance corrections is es- 
timated from the ratio of simulated acceptance distribu- 
tion and its fit function, which is actually used for correc- 
tions. Since the geometry is well understood and a single 
map to exclude malfunctioning areas of the detector has 
been used for the entire data set, this uncertainty is less 
than 1% for all centralities. 

There are two sources of 7r°s not coming from the ver- 
tex (off- vertex 7r°): those produced by hadrons inter- 
acting with detector material (instrumental background) 
and feed-down products from weak decay of higher mass 
hadrons (physics background). Based upon simulations, 
both types of background were found to be negligible at 
less than 1% for px greater than 2.0GeV/c, with the 
exception of 7r°s from K® decay which contribute about 
3% to the 7T° yield for px greater than 1 GeV/c, and have 
been subtracted from the data. The uncertainty due to 
this effect is conservatively estimated as 1.5% and is of 
type C. 



E. 7?aa(A0,pt) 

Similar to the previous analysis [l5[ the Raa{A4>,px) 
measurement uses both the inclusive Raa (jPT ) and the 
quantity v 2 , where v 2 is defined as the second Fourier ex- 
pansion coefficient of the single inclusive azimuthal dis- 
tribution 



dN 



N 



dAcj) 2vr 



(1 + 2v 2 cos(2A</>)) 



(2) 



and Ac/) = "J — (f>. This assumes that the second Fourier 
coefficient is dominant in this expansion. The azimuthal 
anisotropy 1*2 has been published in [16]. 

The 7T° yield is subdivided into six evenly-spaced az- 
imuthal bins in A<f> from to 7r/2 on an event-by-event 
basis using the measured event plane (see Sec. lIIB| ). From 
the inclusive Raa the A(/>-dependent Raa can be con- 
structed as 



Raa(A^, Pt ) = F(A<fi h px) ■ Raa(Pt), (3) 



where 



F(A&,pr) 



N(A(k,PT) 

n 1 

^NiA^px) 



(4) 



and the summation runs over the n — 6 azimuthal bins. 

Because of finite event plane resolution, 
F(A(f>i 7 px) mcas , as calculated from the raw yields, 
needs to be corrected. An approximate unfolding can be 
done by using the raw w™ 1 " an d the resolution-corrected 



y cor 



F(Afr, PT ) = F(A4> uPT ) 



1 



2vf" cos(2A0) 



2wf w cos(2A0) 



(5) 
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The difference between the raw and the corrected «2 is 
given by 



(6) 



(C0S[2(* * - * S )]> 



The denominator is shown in Figure [5J Figure U shows 
the F(Acj),p T ) at 7 < p T < 8GeV/c for centrality 20- 
30%. 
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FIG. 4: (Color online) The corrected ratio, F(A(f>,pT), as a 
function of azimuthal angle at centrality 20-30%. 
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FIG. 6: (Color online) Power law fit parameters (as tabulated 
in Table |nl]) to the 7r° spectra in the 7-20GeV/c pr range 
as a function of centrality, expressed in terms of iVp ar t and 
for p+p (first point). Shown are (a) amplitudes, (b) powers. 
Note that the open points (N pa , It = 352) correspond to 0- 
5 % centrality and partially overlap with the adjacent points 
(iV part =325, 0-10%). 
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FIG. 5: (Color online) Invariant yield of 7r° as a function of 
Pt for each 10% centrality class, 0-5% centrality and mini- 
mum bias. The pr-scale starts at 4 GeV/c. Error bars are the 
sum of statistical and type A systematic uncertainties, boxes 
are the sum of type B and C systematic uncertainties. 



TABLE III: Fit parameters of the power law fit f(pr) = Ap^" 
to the invariant yield (7 < pr < 20GeV/c range) in various 
centrality Au+Au collisions and the p+p cross section [2|| . 



System 


A 


n 


X 2 /NDF 


Au+Au 0-5 % 


no o+ii-b/ 


7.58+0.07 


7.36/9 


Au+Au 0-10 % 


26.3+^ 


7.66+0.05 


5.43/9 


Au+Au 10-20 % 


32.1115 


7.81+0.05 


1.38/9 


Au+Au 20-30 % 


25.6±ag 


7 Q-i+0.06 
'■ 81 -0.05 


14.2/9 


Au+Au 30-40 % 


24.91^ 


7.96+0.06 


11.3/9 


Au+Au 40-50 % 


20.0t 3 3 ;l 


8.02+0.08 


7.50/9 


Au+Au 50-60 % 


15.0+ij-g 


8.09+0.10 


5.56/9 


Au+Au 60-70 % 


5.04±1;£ 


7.92+0.13 


12.6/9 


Au+Au 70-80 % 


ft QO+3.12 


a oq+0.19 

°-°°-o.is 


6.48/8 


Au+Au 80-93 % 




8 7Q+ - 31 


8.14/8 


Au+Au 0-93 % 


16 4 +0 - 93 


7.86+0.02 


11.2/9 


p+p (a) 


16 7+ 173 

lu - ' -1.55 


8.14+0.05 


15.9/9 



Figure [5] shows the ir° invariant yield in Au+Au colli- 
sions for all centralities, and for minimum bias data. As 
with earlier published 7r° results [17] . in this px range 
all distributions are well described by a single power law 
function (f(pr) — Ap7p n ). The fit method employed here 
takes both statistical and systematic uncertainties into 
account, following the one established in previous pub- 
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lications [FjJ, |29|, |30|. The obtained fit parameters are 
listed in Table UTT1 for all Au+Au centrality classes, as well 
as for p+p measured in 2005 [HI]. In the more periph- 
eral collisions the Au+Au and p+p powers are consistent, 
but in central collisions the Au+Au powers are slightly 
smaller, which is also reflected in the behavior of the nu- 
clear modification factor (see Sec. IIII Cj) . In Figure the 
amplitudes and powers from Table Mil are shown. 



B. The production ratio 17/71-° 
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FIG. 7: The 77 to ty° ratio as a function of pr for minimum bias 
Au+Au collisions in 2004 and 2007 (this analysis) data sets, 
dAu collisions (2003), and pythia 6.131 [3l[| . Boxes are pr- 
correlated systematic uncertainties (type B), the small band 
at one is the global uncertainty (type C). 
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FIG. 8: The centrality dependence for the ratio of 77 to n° as 
a function of pt and the expectation from pythia 6.131 [3l| . 
Boxes are pT-correlated systematic uncertainties (type B), the 
small band at one is the global uncertainty (type C). 

Combining the current high statistics ir° results with 
the published 77 meson spectra from the same (2007) data 



FIG. 9: (Color online) One and two sigma standard deviation 
X 2 contours of the linear fit to the minimum bias 77/71° ratio. 



set [301 ] provides new 77/71° ratios with much smaller un- 
certainties than those published previously (32|. In Fig- 
ure [7] the measured 77/71° ratios from minimum bias col- 
lisions are compared for various data sets and colliding 
systems. Although the uncertainties vary, the new ratios 
are consistent with previously published ones [32| and 
are also consistent with the overlaid pythia-6.131 p+p 
calculation. 




0.44 0.46 0.48 0.5 

n/jt (5<p <18(GeV/c)) 



FIG. 10: (Color online) Statistical analysis of the constant fit 
to the minimum bias 77/71° ratio following the method in [29| ]. 

In Figure [8] the 77/71° ratios are shown for various cen- 
tralizes along with the pythia p+p values. A linear 
fit to the minimum bias data gives a constant term of 
0.46+0.05 and a slope of -0.0025+0.0037, with the x 2 
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contours shown in Figure [9l The fit method employed 
here takes both statistical and systematic uncertainties 
into account, following the one established in previous 
publications [l7l . l29l . [301 ] . Since the data are fully consis- 
tent with a zero slope, they were refitted with a constant 
in the 5-18GeV/c pr range, resulting in the final value 
77/71-° = 0.45±0.01(stat)±0.04(syst). Results of the statis- 
tical analysis of the constant fit are shown in Figure [TO] 
Note that the earlier published value [32| for the most 
central Au+Au collisions was 77/71-° =0.40+0.04; the cur- 
rent result is closer to the 77/77° ratios observed in c?Au 
(0.47±0.03) and p+p (0.48±0.03) [H. 

The lack of nuclear effects on this ratio indicate that 
at high pt the fragmentation occurs outside the medium 
and the ratio is governed by vacuum fragmentation [32j . 
This is also supported by a recent global analysis of 77 
fragmentation functions (consider Figure 5 in [33jj and the 
fact that the relevant z range, the fraction of the four- 
momentum of the parton taken by a fragment, in the 
current measurement is about 0.05-0.2). The relevant 
Pt is presumably 5-6GeV/c, below which recombination 
may be a significant hadronization mechanism (see [34l . 
[35l | and [36j]). Also, it should be pointed out, that precise 
knowledge of the absolute value of this ratio is important 
for the background calculations in dielectron and direct 
photon measurements. 



C. Nuclear modification factor ((^-integrated) 

The reference yield of 7r° in P+P collisions has been 
obtained from data taken in 2005 28] . Instead of using a 
fit to the p+p data, Raa has been calculated by dividing 
the Au+Au yields point-by-point by the TAB-scaled p+p 
cross section. Figure [TT] shows Raa for 7r°s as a function 
of pt for six representative centrality classes with the 
new results overlaid on the previously published ones [ijj ■ 
The analysis presented here spans the range pt = 5- 
20GeV/c in several centrality classes. Gray bands show 
the global systematic uncertainties and are of type C, 
which are the quadratic sum of uncertainties of A co n, 
p+p normalization, and off- vertex 7r° contribution shown 
in Table [TTJ The results agree well in the overlapping px 
region with the published Raa data [TtJ ■ 

In Figure [TJ the current RHIC ^/J^ = 200 GeV 
Au+Au 7r° Raa data are compared to the charged 
hadron Raa observed in ^/s^ — 2.76 TeV Pb+Pb col- 
lisions at the LHC (ALICE experiment) [§7|. For both 
centralities and over the entire pr range of 5-20GeV/c, 
the two data sets appear to be similar. This is remarkable 
given the 14-fold increase of colliding energy, resulting in 
an approximately factor of two increase in the parton 
density at the LHC [38J. The expected increase in the 
parton density is corroborated by the factor of 2.2 in- 
crease in dN c ^/d?7 reported by ALICE (39|. 

However, there are two important caveats. Prelim- 
inary results from the same experiment on 7r°s, mea- 
sured via photon conversions up to lOGeV/c, show an 



Raa that is somewhat lower in central collisions than for 
charged hadrons [4(|. In [38[ the authors assert that the 
similarity of Raa at RHIC and LHC may be coincidental. 
In any case, it does not mean that the RHIC and LHC 
data show the same average parton energy loss (e) (see 
Sec. lIIID]) . since the spectra are much harder (the power 
n = 6) at the LHC. The power is obtained by fitting the 
ALICE charged hadron data [37j. 

The fact that at y/s NN = 200 GeV in central colli- 
sions Raa reaches its minimum around 5GeV/c trans- 
verse momentum was first observed in [2|. higher px 
Raa appeared to be approximately constant, although 
the data did not unambiguously exclude a slow rise with 
Pt [H, H3]. On the other hand, all models that repro- 
duce the large suppression observed at pr of 6-10 GeV/ c 
predict a slow rise of Raa as the transverse momentum 
increases [H, |4l|. The current, higher precision data are 
used to reassess the pT-dependence of 7r° suppression in 
the RHIC regime. 

In Figure 1131 a sample linear fit to the pr-dependence 
of Raa in the most central Au+Au data is shown. Fig- 
ure HH shows the 1 and 2er contour lines of the fitted 
slope and intercept for three centralities. The fit method 
employed here takes both statistical and systematic un- 
certainties into account, following the one established in 
previous publications [T3, HH, H(| • In contrast to Figure 
9 in [29( where the slope was consistent with zero within 
lcr due to the large uncertainties, the slope here is signif- 
icantly different from zero, not only in the most central, 
but in 20-30% centrality collisions as well. 

In FigurefTBIthe fitted slopes (a) and the Raa from the 
fits (b) at 7GeV/c and 20GeV/c are shown for all cen- 
tralities, expressed in terms of A part . At and above -/V par t 
=167 (20-30% centrality) the slopes are significantly dif- 
ferent from zero. 



D. Phenomenological energy loss 

The average fractional momentum loss (Si OS s) of high 
Pt hadrons has been of interest since it may reflect the 
average fractional energy loss of the initial parton. S\ oss is 
defined as 5pt/pt, where SpT is the difference of the mo- 
mentum in p+p collisions (pt,pp) and that in Au+Au col- 
lisions (pt.AuAu), and the pt in the denominator is PT,pp- 



In the previous publication [14j . the assumption was 
made that both Au+Au and p+p spectra are comparable 
in shape and Raa vs pt is flat or slowly varying, since 
the data sample size was not large enough to directly cal- 
culate the SpT- With these assumptions, the suppression 
of high pt hadrons could be phenomenologically inter- 
preted as a fractional momentum loss SpT /pt by fitting 
Au+Au spectra with, /(pt) — Ax (pt(1 + ^Pt/pt))^ 71 , 
where A and n were obtained from by fitting a power-law 
function to TAA-scaled p+p cross section [14j . 

With larger statistics p+p and Au+Au data collected, 
it is possible to directly calculate S\ oss without any as- 
sumptions. The calculation method is schematically de- 
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FIG. 11: (Color online) The nuclear modification factor Raa of 7r° as a function of for various 10%-wide centrality classes. 
(Red) closed circles are the results from the current analysis, while (black) open circles are the data published in [13]. Bands 
around 1 indicate global systematic uncertainties and are of type C. The p+p reference is from the 2005 PHENIX data [2^ ]. 



picted in Fig. [16] First, the p+p data is scaled by Taa 
corresponding to the centrality selection of the Au+Au 
data. Next, the p+p point closest in yield to a Au+Au 
point of interest is shifted along the p+p fit function until 
the yield of the p+p and Au+Au points match. The Spr 
is calculated as pr.pp — Pt.AuAu- For obtaining Si oss , the 
5pr is divided by the PT,pp- The uncertainty of the S\ oss 
is calculated by inversely converting the quadratic sum of 
the uncertainties on the yields of Au+Au and p+p points, 
by the p+p fit function. Statistical and type B systematic 
uncertainties are individually calculated in the same way. 
Therefore, the px dependence of systematic uncertainties 
are propagated to the S\ oss values. 

The results are shown in Figure [T7] for minimum bias 
collisions and three different centralities. The uncer- 
tainty coming from Taa, which is of type C, changes 
with centrality selection as listed on the plot. The p+p 
normalization error of 9.7% is not shown here because it 
moves all the points independent of pr or centrality. Be- 
cause S\ oss is plotted as a function of px in p+p collisions, 



the pt points in successive centrality bins in Au+Au are 
shifted as the momentum loss of hadrons varies. An in- 
teresting feature of the central collision data is that while 
5pt/pt is constant up to at least 10 GeV/c, at higher pr 
it slowly decreases, consistent with the slow rise of Raa- 
If one assumes that the fragmentation function of the 
parton after energy loss is unchanged, the fractional mo- 
mentum loss can be interpreted as the average fractional 
energy loss (e) = (AE/E) of the initial parton. This (e) 
can then be compared to the trends predicted in [38[ . In 
this particular model (see Figure 4 in [38]), the collisional 
energy loss appears to be somewhat overestimated, par- 
ticularly below lOGeV/c, but at higher px the observed 
trend in Spt/pt is reproduced quite well. 

In Fig.[T2J it is shown that the Raa in the same central- 
ity at RHIC and LHC show very similar pr dependence 
even though the collision systems and center-of-mass en- 
ergies are vastly different. In Figure 1181 a comparison 
of Slogs is shown. Note that the S\ oss obtained from 
the ALICE charged hadron measurement is ^25% higher 
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FIG. 12: (Color online) Comparison of the 7r° -Raa from this 
measurement and the charged hadron -Raa in Pb+Pb colli- 
sions at ^/s NN = 2.76 TeV from the ALICE experiment [37j 
at LHC. The central and peripheral classes are (a) 0-5% and 
(b) 70-80%. 
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FIG. 13: (Color online) Linear fit to the pr-dependence of 
Raa in the 7-20GeV/c pr range in the most central (0- 
5%) Au+Au collisions. Both statistical (bars) and systematic 
(boxes) uncertainties are considered in the fit. 
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FIG. 14: (Color online) One and two sigma standard devia- 
tion x 2 contours of the linear fit (cf. Figure I13[) to the pr- 
dependence of -Raa for three different centralities. 
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FIG. 15: (Color online) (a) Slopes of the linear fits to it 
Raa vs pr in Js NH =200 GeV Au+Au collisions (as shown 
in Figure I13p and the fitting uncertainties. The fits are in 
the 7-20GeV/c pr range. The horizontal axis is centrality, 
expressed in terms of iVpart- (b) Values of -Raa calculated 
from the fits at pr ~ 7GeV/c and pr = 20GeV/c, also as 
a function of centrality. Note that the open points (7Vp art = 
352) correspond to 0-5 % centrality and partially overlap with 
the adjacent points (N paTt = 325, 0-10%). 



than that from the PHENIX tt° measurement. This is 
reasonable considering the fact that the powers (n) in 
the power-law fit to the px spectra are different between 
the two systems; the power of the PHENIX p+p tt°s at 
= 200GeV/c is about 8, while that of the ALICE 
p+p charged hadrons is about 6. 



E. Model calculations, transport coefficient 

In this section, Raa is compared to four different par- 
ton energy loss models, following the method described 
in [36]. All four models are incorporated into the same 
three-dimensional relativistic hydrodynamic calculation 
with an initial thermalization time tq= 0.6fm/c and de- 
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FIG. 16: (Color online) Method of calculating average frac- 
tional momentum loss (S\ OB s = Spr/pr)- Figure is for illustra- 
tion only, and errors are not shown. In the order of procedure: 
(1) Scale the p+p data by Taa corresponding to centrality se- 
lection of Au+Au data, (2) shift the p+p points closest to 
Au+Au in yield, and (3) calculate momentum difference of 
p+p and Au+Au points. 
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FIG. 17: (Color online) Average fractional momentum loss, 
as defined in the text, between various centrality Au+Au and 
TAA-scaled p+p collisions. The horizontal axis is the pr in the 
p+p collision. Note that for clarity the minimum bias data 
are shifted up by 0.15. <5(global) stands for the uncertainty 
coming from the uncertainties of Taa- The overall normal- 
ization error from the p+p measurement is 1.3%, and is not 
shown here. 
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FIG. 18: (Color online) Comparison of average fractional mo- 
mentum loss, as defined in the text, between the x /s~ A 
200 GeV Au+Au collisions (ty° , current paper) and ./s ? 
2.76 TeV Pb+Pb collisions (ALICE, charged hadrons [37 
The centrality selections are the same. <5(global) stands for the 
uncertainty coming from the uncertainties of Taa- The over- 
all normalization error from the p+p measurement is 1.3% 
for Au+Au data, and is not shown here. 



scribe the observed elliptic flow, pseudorapidity distribu- 
tions, and particle spectra at low px- The Arnold-Moore- 
Yaffe formalism (AMY [9], [42| ) incorporates radiative and 
collisional energy loss processes in an extended medium 
in equilibrium at high temperature, i.e. small coupling 

2 

constant g, where as = f^r- In this approximation, a hi- 
erarchy of scales of successively higher powers of the cou- 
pling constant can be identified, and it becomes possible 
to construct an effective theory of soft modes by sum- 
ming contributions from hard loops into effective propa- 
gators and vertices. The higher- twist approach (HT |l0|) 
is based on the medium-enhanced higher-twist correc- 
tions to the total cross section in deep inelastic scattering 
(DIS) off large nuclei (43|. HT incorporates only radia- 
tive corrections, but it can directly calculate the medium- 
modified fragmentation function. The Armesto-Salgado- 
Wiedemann approach (ASW [ll|), which is equivalent 
to the well known BDMPS-Z approach dHHH, includes 
only radiative processes in a medium where the mean free 
path of the parton is much larger than the color-screening 
length. 

The crucial parameter in all these models is the trans- 
port coefficient q defined as 



q= ^-[GeVVfin] 

A 



(7) 



where /x 2 is the average squared transverse momentum 
transferred from the medium to the parton per collision 
and A is the mean free path of the partons. In AMY q 
is directly related to the temperature, while in HT it is 
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FIG. 19: (Color online) (a) The nuclear modification factor 
of 7T° as a function of pr at centrality 0-5 %. Red solid, green 
dashed and blue dotted-dashed curves are the expectations 
of AMY HIH, HT and ASW El models, respectively, 
(b) The nuclear modification factor of tv° as a function of pr 
at centrality 20-30 %. The theoretical curves in both panels 
are obtained from [36l ]. The gray boxes around 1 show global 
uncertainties and are of type C. 



related to the local entropy density s (oc T 3 ) and in ASW 
it is related to the energy density e. 

In Figure [lj|l the measured nuclear modification fac- 
tor is compared at two centralities with calculations using 
the energy loss models described above, incorporated into 
the same hydrodynamic evolution [36j . In these models, 
the value of q is fixed such as to reproduce the mea- 
sured i?AA in 0-5% centrality collisions. (See (3(| for the 
definitions of the parameters chg, and K, which can be 
converted to <fo.) The values of rfo for gluons (defined as 
the value of q at r = 0.6 fm/c required to describe -Raa) 
differ by a factor of five: <fo is 4.1, 4.3 and 18.5 GeV 2 /fm 
in AMY, HT and ASW, respectively. The HT formal- 
ism was originally developed for deep inelastic scattering 
off a large nucleus, and hence it has become customary 
to quote the value of tfo for a quark [46], and gives the 
value <fo=1.9 GeV 2 /fm as seen in Fig. Despite the 
large differences in the values of q, all models describe 
both the pr- dependence and the centrality dependence 
of -Raa quite well. Additional experimental constraints 
are needed to differentiate between the models, for in- 
stance, restricting the average path length (L) the parton 



traverses in the medium, which can be achieved not only 
in two-particle correlation measurements [l3[ but also by 
studying Raa{A<P) of single particles. 



F. Nuclear modification factor vs event plane 

The overlap region of the colliding nuclei is not az- 
imuthally isotropic, and neither is the medium that is 
formed in the collision. To first approximation (homo- 
geneous density distribution of nucleons) the overlap re- 
gion is elliptical, with the short axis being in the reac- 
tion plane. As a consequence, the average path length 
the hard scattered parton traverses in the medium, los- 
ing energy in the process, varies with the azimuthal an- 
gle A<p, defined experimentally as the relative azimuthal 
angle between the emerging hadron and the measured 
event plane. Measuring Raa as a function of Aip pro- 
vides additional constraints on the average in-medium 
path length [T3 - [T6j . therefore, a more stringent test of 
energy loss models than the ^-integrated Raa alone. 

In Figure [20l the differential nuclear modification fac- 
tor Raa{A<P) is shown for six bins in azimuth and six 
centralities. The participant eccentricities in Table U in- 
dicate the difference between the two extremes, in-plane 
and out-of-plane. In the most central collisions [panel (a)] 
the average path lengths in-plane and out-of-plane are 
almost identical, therefore, the Raa(A<P) curves almost 
completely overlap. As one moves to more peripheral 
collisions, the eccentricity of the overlap region increases 
and the six curves start to split up showing the expected 
ordering: suppression is always largest out-of-plane and 
smallest in-plane. A simple calculation using the partic- 
ipant eccentricity (see Table [I]) shows that the in-plane 
path length changes from 6.1 fm to 3.4 fm when 0-10% 
and 50-60% centralities are compared, while the out-of- 
plane path length changes from 6.7 fm to 5.9 fm between 
the same two centralities. As a consequence, the out-of- 
plane Raa(A4>) changes much less with centrality than 
the in-plane Raa(AiP). All these observations are in full 
agreement with the findings in fl5j ] - 

Figure [21] shows the evolution of Raa(A<P) with cen- 
trality in-plane and out-of-plane at (a) moderate trans- 
verse momenta (averaged in the 6-10GeV/c px region) 
and (b) averaged over all available pr above lOGeV/c. 
As expected, the difference between in-plane and out-of- 
plane suppression increases with eccentricity (decreasing 
-Wpart)) and the actual values converge toward each other 
as the centrality increases. 

Figure [52] shows the comparisons of the models to the 
measured in-plane and out-of-plane -Raa as a function 
of pr for centrality 20-30%. The brackets and bars on 
the data are the statistical uncertainties of the in-plane 
and out-of-plane Raa- The gray band around 1 corre- 
sponds to the systematic uncertainty of the average tt° 
Raa, while the small box at the right end of the line is the 
uncertainty on Taa and is of type C. The bands on the 
data points are the systematic uncertainty of the dN/d(p 
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FIG. 20: (Color online) _Raa(A</!>) as a function of pr for the first six, 10%-wide centrality classes. Each of the six curves 
represent a 15°-wide bin in azimuth, starting from = 0° (in-plane) up to <j> = 90° (out of plane). The gray band around 1 is 
the systematic uncertainty of the normalizing (^-integrated -Raa. The blue boxes around 1 show global uncertainties and are of 
type C. 



including the uncertainty from the event-plane resolu- 
tion. (Red) closed circles and (blue) closed squares are 
the in-plane and out-of-plane Raa, respectively. Panel 
(a) shows the data overlaid with the AMY calculation. 
While the out-of-plane data are well described, the in- 
plane data are not, meaning that the path length de- 
pendence is too weak in this model. The comparison 
with HT on panel (b) shows that this model fails to de- 
scribe both the general trend and the in-plane vs. out- 
of-plane differences. The ASW formalism [panel (c)] de- 
scribes the out-of-plane suppression as well as AMY and 
shows a somewhat larger in-plane vs. out-of-plane dif- 
ference, but still is inconsistent with the data. It should 
be noted that in these three models the energy loss is 
proportional to L 2 , where L is the path length in the 
medium, and the quadratic dependence is characteris- 
tic when radiative energy loss is the dominant mecha- 
nism. Finally the data are compared to a model that 
invokes strong coupling in the medium via an AdS/CFT- 
inspired model. The ASW-AdS/CFT formalism (l2| in- 
corporates the ASW treatment of hard processes, but 
for the soft processes assumes strong coupling. Such a 
hybrid procedure was first suggested in j47|. The vir- 
tual gluons radiated into the medium are governed by 



pQCD, but the interactions of those virtual gluons with 
the medium to bring them on shell is done by assuming 
that the transverse-momentum-squared is proportional 
to L 2 as given by an AdS/CFT calculation [l3.l48j . This 
is in contrast to the the weak-coupling expression for the 
transverse-momentum-squared, qL. This results in an 
energy loss proportional to L 3 instead of L 2 as in the case 
of the pQCD-based models. As can be seen in Figure [22] 
panel (d), the ASW-AdS/CFT model describes both the 
general shape and the absolute difference of the in-plane 
and out-of-plane data well. The observation that mod- 
els with path length dependence of energy loss stronger 
than L 2 are in better agreement with the measurements 
is consistent with the findings in 16]. 



IV. SUMMARY AND CONCLUSIONS 

In summary, the large data set presented in this paper 
made possible a measurement of the 7r° invariant yield 
in ^/s^" = 200 GcV Au+Au collisions up to 20GeV/c 
transverse momentum. This has led to a precision mea- 
surement of the i]/ir ratio in Au+Au collisions, which is 
constant as a function of both centrality and px, r)/ir a — 
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FIG. 21: (Color online) Centrality dependence (expressed 
in terms of Af part ) of the ir° _Raa(A</>) in-plane (red, upper 
points) and out-of-plane (blue, lower points), averaged (a) in 
the 6-10GeV/c transverse momentum region and (b) above 
lOGeV/c. Boxes are systematic uncertainties on 7?aa(A0). 



0.45±0.01(stat)±0.04(syst) and consistent with the val- 
ues observed in dAu and p+p. The large observed ir° 
suppression is fully consistent with earlier findings, and 
a slow but significant rise of Raa vs pt with a slope 
of O.OlOei^KGeV/cT 1 for central collisions is now 
observed for the first time at RHIC energies. This has 
been an expectation of all pQCD-based parton energy 
loss models. The large data set has also made pos- 
sible the calculation of a phenomenological AE/E en- 
ergy loss. The differential i?AA(A</>), testing the path- 
length dependence of energy loss, is measured up to px 
of 20GeV/c and is compared to various energy loss cal- 
culations. While all models considered describe the <j>- 
integrated Raa adequately, the pQCD based calculations 
where the energy loss depends on the path length as L 2 
fail to describe the differential Raa(A(/>). The data re- 
quire an energy loss with a power greater than 2, as 
given by models in which the soft interactions with the 
medium are strongly coupled. These findings are consis- 
tent with the conclusions of .16] in which data on elliptic 
flow of high pt- (> 6 GeV) 7r°s is shown to be incon- 
sistent with pQCD-based models. To explore the strong 
coupling regime, a comparison was made to the same 



ASW-AdS/CFT model used in this work, as well as to 
a phenomenological model [49| in which the energy loss 
was proportional to L 3 , both of which were able to fit the 
data. Both the current measurement and [l6| explore a 
region of high py where the mechanism leading to an az- 
imuthally anisotropic yield is parton energy loss rather 
than hydrodynamical flow. It is increasingly difficult for 
purely pQCD-based models to explain these results and 
one is led to the tentative conclusion that strong cou- 
pling plays an important role in parton energy loss in 
the medium. At present, the best method to do the rele- 
vant calculations is in an AdS / CFT framework. Similar 
conclusions are reached when one looks at the behav- 
ior of heavy quarks [50j , where higher quality data will 
soon be available. Recent preliminary results on the sup- 
pression pattern seen at the LHC for py > 6 GeV/c are 
strikingly similar to those seen at RHIC. In this paper, 
a phenomenological calculation of fractional energy loss 
is given, which indicates that the energy loss at LHC 
(ALICE data) is about 25% higher than at RHIC, and 
that the loss falls slightly with energy. The dependence 
of these observables on momentum and center-of-mass 
energy (presumably on energy density) will be a crucial 
factor in untangling the underlying mechanisms of parton 
energy loss. Recently, experiments at the LHC have be- 
gun to examine the behavior of fully reconstructed jets, 
which should give more easily interpretable information 
on this phenomenon. Future work at both the LHC and 
at RHIC should bring data on path length dependence of 
fully reconstructed jets, jet widths, and heavy-quark jets 
which will add a wealth of new information. In addition, 
a more complete understanding of the initial state is also 
needed both for the initial configuration of hydrodynam- 
ical models, and as a calibration of the hard probes that 
are used in these measurements. 
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FIG. 22: (Color online) The data points are 7?aa(A0) in 20-30% centrality as a function of pr for in-plane (red, upper points) 
and out-of-plane ir° (blue, lower points), compared to four model calculations (see text for description and references), (a) 
pQCD-based AMY [a El, (b) HT [nj, (c) pQCD-based ASW [TH, (d) ASW using AdS/CFT correspondence [H. The curves 
on panels (a)-(c) are taken from[36(|. The dotted and solid lines are the in-plane and out-of-plane predictions, respectively. The 
definition of the gray bands and blue boxes is same as in Fig. I20I 
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APPENDIX 

Tables IIVIIVII give values for the invariant Yields for 
neutral pions, as shown in Fig. [SJ Tables IVIIHIXI give 
values of Raa for neutral pions, as shown in Fig. ITT1 



TABLE IV: Invariant yields of neutral pions as a function 
of pr at \y\ < 0.35 in Au+Au collisions at - v /ijviv=200 GeV 
for the very- most central 0-5% centrality. Syst.(B) refers to 
type-B systematic errors. See Fig. [5] 
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TABLE V: Invariant yields of neutral pions as a function of 
Pt at \y\ < 0.35 in Au+Au collisions at - v /sjvjv=200 GeV 
for 0-10%, 10-20%, 20-30%, 30-40%, and 40-50% central- 
ity ranges. Syst.(B) refers to type-B systematic errors. See 
Fig-El 



TABLE VI: Invariant yields of neutral pions as a function 
of pt at \y\ < 0.35 in Au+Au collisions at v / 5jvjv=200 GeV 
for 50-60%, 60-70%, 70-80%, 80-93%, and 0-93% central- 
ity ranges. Syst.(B) refers to type-B systematic errors. See 
Fig. E] 
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TABLE VII: Nuclear modification factors, Raa for neutral 
pions as a function of pt at \y\ < 0.35 in Au+Au collisions at 
y/sj^=200 GeV for 0-10%, 10-20%, 20-30%, 30-40%, and 
40-50% centrality ranges. Syst.(B) refers to type-B system- 
atic errors. The global systematic uncertainties (type C) are 
p+p normalization (9.7%) and Off-vertex (1.5%). See Fig.fTTI 



TABLE VIII: Nuclear modification factors, -Raa for neutral 
pions as a function of pr at |yj < 0.35 in Au+Au collisions at 
■ s /sM=20Q GeV for 50-60%, 60-70%, 70-80%, 80-93%, and 
0-93% centrality ranges. Syst.(B) refers to type-B systematic 
errors. The global systematic uncertainties (type C) are p+p 
normalization (9.7%) and Off-vertex (1.5%). See Fig. [Til 
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33 


0.85 


36 


19 





.3240 


0.13 


41 


0.14 


44 






19 





8040 


0.82 


100 


0.35 


44 


5.25 





.3528 


0.0024 


0.69 


0.045 


13 


70 


-80% 


5.25 





.8701 


0.014 


1.6 


0.11 


13 


5.75 





.3377 


0.0031 


0.92 


0.043 


13 






5.75 





8731 


0.019 


2.2 


0.11 


13 


6.25 





.3412 


0.0041 


1.2 


0.044 


13 






6.25 





8736 


0.026 


3.0 


0.11 


13 


6.75 





.3382 


0.0053 


1.6 


0.044 


13 






6.75 





8608 


0.034 


3.9 


0.11 


13 


7.25 





.3347 


0.0065 


2.0 


0.043 


13 






7.25 





9298 


0.045 


4.9 


0.12 


13 


7.75 





.3444 


0.0083 


2.4 


0.045 


13 






7.75 





8904 


0.059 


6.6 


0.12 


13 


8.25 





.3511 


0.010 


3.0 


0.043 


12 






8.25 





8823 


0.069 


7.8 


0.11 


12 


8.75 





.3268 


0.012 


3.6 


0.040 


12 






8.75 





9520 


0.083 


8.7 


0.12 


12 


9.25 





.3425 


0.015 


4.4 


0.042 


12 






9.25 





8898 


0.11 


12 


0.11 


12 


9.75 





.4006 


0.021 


5.1 


0.051 


13 






9.75 


1 


.024 


0.13 


13 


0.13 


13 


11 





.3779 


0.014 


3.7 


0.057 


15 






11 


0.7552 


0.079 


10 


0.11 


15 


13 





.4085 


0.028 


6.9 


0.085 


21 






13 


1.069 


0.18 


17 


0.22 


21 


1") 





.5216 


0.065 


12 


0.15 


28 






15 


0.6448 


0.29 


45 


0.18 


28 


17 





.3713 


0.095 


26 


0.13 


30 






17 


1.023 


0.74 


72 


0.37 


36 


19 





.4550 


0.20 


43 


0.20 


44 
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5.25 





.4408 


0.0032 


0.72 


0.056 


13 


80 


-93% 


5.25 





8298 


0.021 


2.5 


0.11 


13 


5.75 





.4234 


0.0041 


0.98 


0.054 


13 






5.75 





8385 


0.028 


3.4 


0.11 


13 


6.25 





.4394 


0.0056 


1.3 


0.056 


13 






6.25 





.8537 


0.039 


4.0 


0.11 


13 


6.75 





.4341 


0.0072 


1.7 


0.056 


13 






6.75 





9070 


0.051 


5.6 


0.12 


13 


7.25 





.4186 


0.0089 


2.1 


0.054 


13 






7.25 





9464 


0.067 


7.1 


0.12 


13 


7.75 





.4398 


0.011 


2.6 


0.057 


13 






7.75 





8566 


0.081 


9.5 


0.11 


13 


8.25 





.4397 


0.014 


3.2 


0.053 


12 






8.25 





9038 


0.11 


12 


0.11 


12 


8.75 





.4313 


0.017 


3.9 


0.053 


12 






8.75 





6904 


0.13 


19 


0.084 


12 


9.25 





.4193 


0.020 


4.9 


0.052 


12 






9.25 





7566 


0.14 


18 


0.093 


12 


9.75 





.4291 


0.025 


5.9 


0.054 


13 






9.75 





.6419 


0.15 


24 


0.081 


13 


11 





.4526 


0.019 


4.1 


0.068 


15 






11 





9396 


0.13 


14 


0.14 


15 


13 





.4762 


0.036 


7.6 


0.099 


21 






13 





6637 


0.21 


32 


0.14 


21 


15 





.5969 


0.092 


15 


0.17 


28 






15 





2963 


0.30 


100 


0.083 


28 


17 





.5002 


0.14 


28 


0.18 


36 






17 


1 


.176 


1.2 


100 


0.42 


36 


19 





.5423 


0.27 


49 


0.24 


44 
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5.25 





.5422 


0.0043 


0.79 


0.069 


13 


0- 


93% 


5.25 





.3105 


0.0014 


0.46 


0.040 


13 


5.75 





.5346 


0.0057 


1.1 


0.069 


13 






5.75 





3002 


0.0019 


0.62 


0.038 


13 


6.25 





.5532 


0.0078 


1.4 


0.071 


13 






6.25 





.3013 


0.0025 


0.82 


0.039 


13 


6.75 





.5531 


0.010 


1.8 


0.071 


13 






6.75 





3025 


0.0032 


1.1 


0.039 


13 


7.25 





.5227 


0.012 


2.4 


0.067 


13 






7.25 





2962 


0.0040 


1.4 


0.038 


13 


7.75 





.5562 


0.016 


2.9 


0.072 


13 






7.75 





.3031 


0.0051 


1.7 


0.039 


13 


8.25 





.5404 


0.020 


3.6 


0.066 


12 






8.25 





3070 


0.0064 


2.1 


0.037 


12 


8.75 





.5457 


0.024 


4.3 


0.067 


12 






8.75 





2973 


0.0074 


2.5 


0.036 


12 


9.25 





.5466 


0.029 


5.3 


0.067 


12 






9.25 





.3123 


0.0096 


3.1 


0.038 


12 


9.75 





.6059 


0.038 


6.3 


0.077 


13 






9.75 





3368 


0.013 


3.8 


0.043 


13 


11 





.5231 


0.025 


4.8 


0.078 


15 






11 





3240 


0.0088 


2.7 


0.048 


15 


13 





.5690 


0.051 


9.0 


0.12 


21 






13 





.3458 


0.018 


5.1 


0.072 


21 


15 





.4362 


0.25 


57 


0.12 


28 






15 





.4370 


0.040 


9.2 


0.12 


28 


17 





.6280 


0.20 


32 


0.22 


30 






17 





.4344 


0.067 


10 


0.16 


30 


19 





.9655 


0.48 


49 


0.42 


44 






19 





3750 


0.10 


27 


0.16 


44 
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TABLE IX: Nuclear modification factors, Raa for neutral pi- 
ons as a function of pr at \y\ < 0.35 in Au+Au collisions at 
V s iViv = 200 GeV for the very most central, 0-5% collisions. 
Syst.(B) refers to type-B systematic errors. The global sys- 
tematic uncertainties (type C) are p+p normalization (9.7%) 
and Off-vertex (1.5%). See Fig. ITU 



Ccn- 


Pt 


-Raa 


Stat. 


Fraction 


Syst.(B) 


Fraction 


trality 






error 


% 


error 


% 


0-5% 


5.25 


0.1753 


0.0018 


1.0 


0.022 


13 




5.75 


0.1753 


0.0022 


1.3 


0.022 


13 




6.25 


0.1756 


0.0028 


1.6 


0.023 


13 




6.75 


0.1823 


0.0036 


2.0 


0.023 


13 




7.25 


0.1749 


0.0043 


2.4 


0.023 


13 




7.75 


0.1815 


0.0053 


2.9 


0.023 


13 




8.25 


0.1894 


0.0067 


3.5 


0.023 


12 




8.75 


0.1766 


0.0076 


4.3 


0.022 


12 




9.25 


0.2080 


0.010 


4.8 


0.026 


12 




9.75 


0.2288 


0.013 


5.7 


0.029 


13 




11 


0.2079 


0.0087 


4.2 


0.031 


15 




13 


0.2455 


0.018 


7.3 


0.051 


21 




15 


0.2982 


0.042 


14 


0.083 


28 




17 


0.3137 


0.076 


24 


0.11 


36 




19 


0.3308 


0.14 


43 


0.14 


44 
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